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Q\ . HST WFPC2 archive, of which 43 are type SO or earlier. In this paper, we present the 

(~| , data set and derive the V — I color distributions. We derive the first four moments of 

*^| the color distributions, as well as a measure for their non-unimodality. The number of 



S 



We have compiled data for the globular cluster systems of 50 galaxies from the 



O I globular clusters in each galaxy ranges from 18 (in NGC 2778) to 781 (NGC 5846). For 

'^ ' those systems having more than 100 clusters, seven of sixteen (44%) show significant 

C^ ■ bimodality. Overall, roughly half of all the systems in our sample show hints of a 

bimodal color distribution. In general, the distributions of the faint galaxies are 
^ . consistent with unimodality, whereas those of the brighter galaxies are not. We also 

find a number of systems with narrow color distributions — with both mean red and 
blue colors — suggesting that systems exist with only metal-rich or only metal-poor 
globular clusters. We discuss their possible origins. 

In comparing the moments of the V — I distributions with various galaxy properties 
for the early-type galaxies, we find the following difference in the correlations between 
the field and cluster galaxy populations: the peak V — I color of the globular cluster 
distribution correlates well with the central velocity dispersion — and hence the Mg2 
index and total luminosity — for galaxies in cluster environments; there exists no such 
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correlation for field galaxies. This difference between cluster and field galaxies possibly 
reflects different formation scenarios for their globular cluster systems. Among the 
explanations for such a correlation, we consider either a larger age spread in the field 
populations or the possibility that cluster galaxies are always affected by significant 
accretion whereas some field galaxies could host pure "in situ" formed populations. 



Subject headings: galaxies: star clusters, galaxies: elliptical and lenticular, cD 



1. Introduction 

In the past decade, globular cluster systems were established as powerful probes for the 
formation and evolution of galaxies (e.g. Harris 1991 and Ashman & Zepf 1998 for reviews). With 
very few exceptions, globular clusters exist in large numbers around all galaxies. They are not 
only easy to detect; their properties are readily influenced by the formation and evolution of their 
host galaxies and, in a larger sense, their environment. Globular cluster formation is also closely 
related to the star-formation history of the galaxies (e.g. Kissler-Patig et al. 1998a). The study 
of globular cluster systems, combined with the predictions of their properties from various galaxy 
formation scenarios, can therefore be used to constrain galaxy formation and evolution. 

In recent years, tens of galaxies have been studied for their globular cluster population 
(e.g. compilation in Ashman &: Zepf 1998). However, no systematic study of the properties of 
globular cluster systems, including more than a dozen hosts, has been carried out. Motivated by 
the need for a homogeneous study of a large number of systems, we have exploited the Archive of 
the Hubble Space Telescope as an ideal starting point for such a study, since it includes a large 
number of images centered on galaxies in different filters and taken in nearly identical conditions, 
leading to homogeneous photometry. 

We begin by studying the properties of globular cluster systems around a large number of 
early-type galaxies, then describe the resulting sample in Sect. 2 and the pipeline reduction in 
Sect. 3. We study the color distributions of the systems (presented in Sect. 4 and 5), followed by 
a discussion of the properties of the color distributions, correlations among these properties, and 
their relation to the properties of the host galaxies (Sect. 6 and 7). Because color denotes a mix 
of age and metallicity of the systems, globular cluster systems have gained considerable interest 
after the discovery of multi-modal distributions in some systems (Zepf &: Ashman 1993). Such 
a system implies that at least in some galaxies, globular clusters formed from several different 
mechanisms and possibly at several epochs. Color distributions alone are thus already a strong 
diagnostic for the history of the globular clusters systems and, by implication, for the history of 
the host galaxies. We discuss the properties of the color distributions of our sample for various 
formation scenarios in Sect. 8. In a companion paper, we use these data to constrain the models 
of Cote et al. (1998). 



2. The Galaxy Sample 

We searched the HST archive for all early-type galaxies that have been observed in two 
colors with adequate exposure times, using only WFPC2 data from a variety of programs, with 
the Nuker Group (PI: Faber and Richstone) providing most of it. 90% of the galaxies where 
observed in F555W (V) and F814W (/); some galaxies were observed through the F547M and the 
F450W filters which we transformed to V (see below). We did not use data taken in other filters. 
Although we may have some some early-type galaxies, we attempted to make a complete study 
from the archive. Furthermore, some of the galaxies below are not early-types; we include them 
only for comparison, not for use in the analysis. 

Table 1 presents the data properties for each galaxy: galaxy name (Col. 1), filter (Col. 2), 
total exposure time (Col. 3), number of individual exposures (Col. 4), and PI name and number 
(Col. 5). Three galaxies have more than one pointings and we include their total exposure times 
from all pointings. There are 50 galaxies in total. Table 2 summarizes the physical properties for 
the 43 galaxies which are SO or earlier. 



3. The Pipelined Reduction 

We reduced all of the data starting from the individual exposure. The CADC provided the 
recalibrated science exposure (see cadcwww.dao.nrc.ca for details of the recalibration) . All data 
were sent through a pipelined reduction. We divide the reduction into two steps: the combination 
of the individual images, and the photometry for the globular clusters. The first step requires 
shifting the images, scaling by exposure times, and combining. We use the HST header parameters 
for the telescope pointing (CRVALl and CRVAL2) to determine the shift. Instead of using the 
telescope pointing, we could have used reference objects in the images to measure the shift; for a 
handful of galaxies, we confirmed that the telescope pointing provides the same shift as when we 
used reference objects. We use a linear shift so as not to smear cosmic-ray affected pixels, then 
scale the images to a common exposure time. Different gains were used between galaxies and even 
among different exposures of the same galaxy. Finally, we transformed all images to gain 7. The 
combining process must be robust, as many of the galaxies have a limited number of individual 
exposures that make it difficult to remove cosmic-rays adequately. We therefore use a biweight 
estimate for the combining. The biweight attempts to estimate the mode of a distribution by 
using an iterative approach whereby deviant points are subsequently de-weighted. The parameters 
for the weighting result from numerical simulations (see Beers et al. 1991 and reference therein for 
a fuller explanation). Unfortunately, some galaxies only have two exposures in a given filter. To 
help decide whether a given pixel is an outlier, we include both filters in the combining process. 
We exclude the flagged pixels when combining in a given filter. 

The second part of the reduction includes the photometry. Most of these data have the 
center of the galaxy in one of the WFPC2 chips (most are in the PC). To measure with greater 



reliability those clusters which have the high-background galaxy under them, we first subtract 
the galaxy light using a "median-filtering" technique; instead of the median, we use a biweight 
estimate (explained above) for the local flux value. We choose the parameters for the filtering — i.e., 
the inner radius equals 8 pixels and the outer radius equals 14 pixels — based on results from 
experiments with no background present. In addition, we make an initial pass to detect the 
globular clusters, then flag those pixels under the clusters to make sure the filtering technique 
does not include them when determining the galaxy component. We subtract the galaxy only in 
the chip that contains the galaxy center. 

We use SExtractor (Bertin & Arnouts 1996) for the photometry, setting the finding parameters 
to 2 connected pixels 2.5 a over the local background as defined in a 32x32 pixels box. The frames 
were further convolved with the PSF kernel given in the WFPC2 handbook for the respective 
chips in order to optimize the finding (see SExtractor manual, Bertin & Arnouts 1996). 

We must correct the resulting magnitudes since the photometry is measured in a 2 pixel 
radius aperture. We correct the resulting magnitudes to a 0.5" radius aperture, using the values 
proposed in Whitmore et al. (1997), and by 0.1 magnitude in order to extrapolate to the total 
light (see Holtzman et al. 1995). Finally, the magnitudes measured in the HST filter set were 
converted into the Johnson-Cousins V and / system by applying the transformations given in 
Holtzman et al. (1995). For the F555W and F814W filters used in the vast majority of our data 
sets, empirical transformations exists. For the F547M filter, we used the theoretical transformation 
into the V filter. Finally, we used the empirical transformation to convert the F450W filter, 
used for the Fornax galaxies, into a Johnson B filter, followed by a linear transformation from 
B — I into V — I derived from the populations synthesis models of Bruzual & Chariot (1993) 
for ages > 5 Gyrs and the full range of available metallicities. The linear transformation follows 
V — I = 0.466 {B — I) + 0.14. The average uncertainty on the magnitudes is 0.06; this uncertainty 
scales with magnitude roughly so that at V=2A we have an uncertainty of 0.1, and at V=2Q an 
uncertainty of 0.01. We go down to around 25.5 in V and 24.5 in /, but we are only able to 
approximate both our limiting magnitude and uncertainty, given that the galaxies have a range of 
exposure times. 

We note that the extrapolations from 2 pixel apertures to total magnitudes, derived for point 
sources, are extremely uncertain when applied to slightly extended objects such as globular clusters 
and can significantly affect the final magnitudes (see Kundu & Whitmore 1998, Puzia et al. 1999). 
However, the error of the extrapolation is similar in the F555W and F814W filters, so that the 
final colors are not affected by systematic errors of less than 0.03 magnitudes. 

Finally, the data were corrected for foreground reddening using the maps of Schlegel et al. 
(1998), and the law E(y - /) = 1.55 -{B-V). 



4. The V — I color distribution of the systems 

SExtractor provides an estimate of the magnitude, magnitude error, ellipticity, FWHM, 
and star/galaxy classification. From these measurements we select globular clusters from the 
objects that have the following criteria: a magnitude error smaller than 0.1; V — I color between 
and 2; ellipticity smaller than 0.5; FWHM between 1 and 4 pixels (on both the PC and WF 
chips); and 'non-star' classification (set to reject remaining cosmic rays and hot pixels, rather 
than to try to exclude point sources). No cut on the magnitude exists since it correlates with the 
magnitude error. The number of true globular clusters excluded by these criteria is typically less 
than 10%, somewhat higher for the most distant galaxies and those with short exposure times. 
The contamination by foreground stars and background galaxies passing our selection criteria is 
roughly constant and estimated from the Hubble Deep Field North (Williams et a/. 1996). The 
percentage of contamination depends on the distance and the richness of the globular cluster 
system considered, and lies around 5% for a typical system with > 100 clusters at a distance of 
10 to 20 Mpc. It can rise up to 30% to 50% for systems at significantly larger distances — that we 
then excluded in the following analysis where mentioned — or for systems with significantly less 
clusters. 

Fig. 1 plots V — I histograms for all of the 50 galaxies in the sample. As a histogram is the 
crudest form of density estimation, we also include the adaptive kernel density estimate as the 
solid line. A histogram is a density estimate that utilizes both discreet and equal-sized windows; 
it is therefore dependent on both the choice of the window size and the location of the window 
edges. More robust techniques include both weights for an individual datapoint's contribution to 
a density estimate, and varying size that preserves the noise properties (i.e., sparse regions require 
larger binning). The adaptive kernel estimate contains these properties, provides a non-parametric 
density estimate on a grid. We refer to Silverman (1986) for a complete discussion of the adaptive 
kernel estimate. It works as follows: we obtain an initial density estimate, then adjust the 
window width at each grid element according to the initial density in such a way that low density 
regions have a large window width. We use the Epanechnikov kernel, an inverted parabola, for 
the density estimate. The result of an adaptive kernel estimate can be highly dependent on the 
initial smoothing, especially for small samples like ours with N < 100; thus, we use least-squares 
cross-validation to choose the initial smoothing (see Silverman 1986, p. 48). This technique 
minimizes an estimate of the integrated square error given by /(/ — /)^, where / is the density 
estimated from the data — dependent on the window width (the smoothing parameter) — and / is 
the underlying function we are trying to estimate. One point is removed from the sample, and we 
estimate the density for a given smoothing from the N — 1 data points. Repeating this N times 
and using the average provides an estimate of the underlying function /. The optimal smoothing 
is then the one that minimizes the estimated integrated square error. 

The dotted lines in Fig. 1 represent the 90% confidence bands obtained through a bootstrap 
estimate, a technique that draws with replacement from the original sample of N points. One 
bootstrap realization contains N points, which are not necessarily unique. For each realization we 



estimate the distribution function and repeat this procedure 1000 times. At each grid element we 
then have a distribution of points, from which we choose the 5% and the 95% values, representing 
the 90% confidence band. These confidence bands reflect only the variance, not the bias, of the 
estimate. The bootstrap procedure includes the measurement uncertainties implicitly because 
they are reflected in the original estimate of the colors. 



4.1. Measuring the Moments of the Distributions 

Moments represent a simple quantification of the distributions to with which we correlate other 
galaxy properties (discussed below). However, due to small number statistics and contamination, 
we desire a more robust estimate for the moments other than the classical estimates. We therefore 
use the biweight estimate of location (1st moment) and scale (2nd moment). Finch's asymmetry 
index (3rd moment), and the Tail Index (4th moment). Beers et al. (1991) describe the biweight 
estimators in detail. The asymmetry index measures skewness by comparing gaps between the 
ordered data points on both the right and left side of the distribution. The gaps are weighted 
according to their rank order, thereby devaluing the contribution from outliers. The Tail Index 
measures the spread of the dataset at the 90% level relative to the spread at 75%, and is normalized 
such that a Gaussian has an index equal one. Bird & Beers (1993) provide an explanation of 
the asymmetry index and tail index, and show that they are superior estimates compared to the 
classical skewness and kurtosis. Table 3 presents these robust moments and their 68% confidence 
band. A bootstrap procedure, the same as described above, provides the confidence bands. The 
galaxies listed in italics in Table 3 are late-type. 



5. Relations between Moments and Galaxy Properties 

In the following sections, we investigate the relations among the properties of the globular 
cluster color distributions and their relation to the properties of the host galaxies. We restrict 
ourselves to a homogeneous sample of globular cluster systems with best determined properties. 
We include SO's and ellipticals, but exclude the late-type galaxies listed in italics in Table 3. We 
further restrict the sample to galaxies with redshift less than 3000 km s~^ (excluding NGC 1700 
and NGC 4881), since at anything smaller than this distance contaminating fore-/background 
objects present no problem. We set a lower limit of the number of globular clusters of 50, for 
which the moments are well defined. Finally, we exclude the systems that were not observed in 
the F555W and F814W filters to avoid systematic effects in the conversion to a homogeneous 
photometric system (excluding NGC 4374, NGC 1399, NGC 1404). 

The restricted sample includes 26 globular cluster systems for which the distributions of 
properties are shown in Figure 2. Figure 2 presents all properties plotted against each other, as 
we discuss below. Examining the range of the V — I moments shows a relatively sharp cut-off in 



the peak color at F — / ~ 1.12. The width or scale appears to be around 0.16 for most systems, 
including an artificial broadening by photometric errors of around 0.06 typically, with only few 
systems having significantly broader distributions. Most of the systems are slightly positively 
skewed — i.e., they have the tendency to cut-off more sharply in the blue and extend somewhat 
towards the red. Finally, the tail index implies that most of the distributions have tail weight 
consistent with a Gaussian distribution. 



5.1. Peak to the 2nd, 3rd, and 4th Moments 

We show the relation between the peak and the scale (width) of the distributions in the first 
plot on the upper line in Figure 2; an enlarged version appears in Figure 3. We first notice that a 
blue cut-off is seen in all distributions, implying a lower metallicity limit for all globular clusters. 
This limitation imposes a lower cut-off in color and forbids a broadening of any distribution 
towards the extreme blue (about V — I < 0.8). Therefore, in order to have a blue mean color, 
systems both must lack those red globular clusters that could not be compensated by extreme blue 
clusters and must be narrow. Interestingly, such systems with very little red clusters are observed 
in our sample, i.e. seem to exist among early-type galaxies. 

On the other hand, red systems show a wide range of scale values, indicating that both systems 
with only red clusters (red peak, small scale) and systems with a wide range of globular cluster 
colors exist. Figure 3 shows that the reddest systems (about V — I > 1.1) have systematically large 
scales. These systems must have a broad range of colors, including blue clusters. For systems with 
V — I > 1.1, a large scale is partly expected due to the fact that the V — I color changes faster 
with metallicity at high metallicities (e.g. Kissler-Patig et al. 1998b, compare with the theoretical 
predictions, e.g. Worthey 1994) and thus broadens any distribution. 

The comparisons of the peak color to either the skewness or the tail indicate no significant 
trends. The exception is expected: if a blue system is strongly skewed, the skewness is towards 
the red color, and vise-versa. However, not only do most systems fail to show a strong skewness, 
the tail weights of the color distributions do not deviate much from a Gaussian and show no 
systematic trends. 



5.2. Properties of the host galaxies 

For the host galaxies we compiled a number of properties to correlate with the properties 
of the globular cluster color distribution. These galaxy properties are given in Table 2, together 
with the source. The galaxy properties include type, total luminosity, total color, central velocity 
dispersion, central Mg2 index, environment density, isophotal shape parameter, and, when 
applicable, effective radius. 



The distribution of the properties for the 26 galaxies appear in the plots in the lower-right 
side of Figure 2. The vast majority are elliptical galaxies. We separate field from cluster galaxies 
at a density parameter of p = 1: i.e., galaxies in groups such as Leo (NGC 3377, NGC 3379, 
NGC 3384) will belong to the field, while galaxies in small or medium sized clusters, such as 
Fornax and Virgo, will belong to the cluster category. We thus have 15 galaxies belonging to the 
field (the green triangles in Figure 2) and 11 belonging to the cluster category (the blue circles). 

While the Mg2 - a relation is well defined for our galaxies, the Faber-Jackson relation shows 
significant scatter, probably due to noise in the distances. Therefore, we use the central velocity 
dispersion as a distance independent measure for the relative galaxy mass, rather than the absolute 
magnitude. We have also made comparison using log(re o"^), which should reflect the galaxy mass 
better than the dispersion alone, and find no better correlations. 



5.3. Correlations with the Whole Sample 

We first consider the whole sample, independent of environment. No significant correlation 
between the peak or scale of the globular cluster systems and any galaxy property could be 
detected, as measured by the Spearman rank-order correlation coefficients. We only note that 
boxy galaxies (a/a4 < 0) have systematically red peaks, and that the peak values scatter most for 
intermediate galaxies when plotted against My, a, or Mg2 of the host galaxies. 

At face value, this result might contradict the known dependence of the mean metallicity of 
globular cluster systems from galaxy luminosity (e.g. Harris 1991, Forbes et al. 1997), in so far 
as brighter galaxies host globular clusters with a higher mean metallicity. Therefore, we could 
have expected a systematically redder peak with increasing My or a of the host galaxy. But 
Ashman & Zepf 1998 noticed that the mean color of the clusters does not correlate with the galaxy 
luminosity for early-type galaxies. In fact, the relation between mean cluster metallicity and host 
galaxy luminosity (cf. Harris 1991) results from the fact that dwarf galaxies possess less metal-rich 
globulars than spirals, which in turn host less metal-rich globulars than giant early-type galaxies. 
Within a given galaxy type, the relation between mean cluster metallicity and galaxy luminosity 
is much less clear, if present at all. 



5.4. Comparing Field and Clusters Galaxies 

The recent models of hierarchical clustering formation for galaxies (Kauffmann 1996, 
Baugh et al. 1996) predict significant differences in the formation and star-formation histories for 
field and cluster galaxies. For this reason, we investigate these two classes separately and split our 
samples of galaxies into field and cluster galaxies, as described above. 

According to a Spearman rank-order correlation coefficients, no correlations are detected 



in the field sample of 15 galaxies. However, for cluster galaxies, the peak values are found to 
correlate positively with the host galaxy velocity dispersions at the > 99% confidence level and, 
accordingly, at > 98% confidence level with their Mg2 index values. Peak against galaxy velocity 
dispersion appear in the seventh plot on the upper-line of Figure 2 for field and cluster galaxies; 
an enlarged version is shown in Figure 3. The same results are true when using the total galaxy 
magnitude — i.e., cluster galaxies show a correlation, whereas the field galaxies do not. 



6. Bimodality 

Tests for bimodality are difficult to quantify. In many of the past globular cluster systems 
studies, authors have used double Gaussian fits to decide whether the system is bimodal. However, 
these results are highly dependent both on the underlying assumption of Gaussianity and on the 
criteria for deciding the number of sub-populations. More general tests exist in the statistics 
literature that impose global conditions on unimodality with no underlying assumptions about 
the sub-populations. Here, we use the Dip statistic (Hartigan & Hartigan 1985). The Dip test 
measures the maximum distance between the empirical distribution function and the best fitting 
unimodal distribution, with the significance level based on simulations. Table 3 lists the value for 
the significance of the Dip statistic, but not the Dip value itself. The same is true in Figure 2 and 
4 where we plot the significance value only. A significance of 0.9 implies a 90% probability that 
the distribution is not unimodal. We do not calculate values below 0.5, as the Dip statistic is only 
defined to discriminate against unimodality. 

Figure 4 presents the DIP probabilities and confidence bands for all of the galaxies. Out of 
43 galaxies for which the dip can be measured, 9 galaxies show clear deviation from unimodality. 
These are NGC 1399, NGC 1404, NGC 1426, NGC 3640, NGC 4472, NGC 4486, NGC 4526, 
NGC 4649, and NGC 5846. Six of these were already known from ground or space studies to show 
several distinct populations (NGC 1399: Forbes et al. 1998, Ostrov et al. 1998, Kissler-Patig et al. 
1997a; NGC 1404: Forbes et al. 1998; NGC 4472: Geisler et al. 1996, Puzia et al. 1999; 
NGC 4486: Whitmore et al. 1995, Elson & Santiago 1996; NGC 4649: Neilsen et al. 1999; 
NGC 5846: Forbes et al. 1997). Using a test to discriminate between single and double Gaussians, 
Neilsen et al. (1999) found significant bi-modal population in four of their eight galaxies studied. 
Using the same datasets, we confirm most of their measurements of bi-modality; we disagree only 
in NGC 4552 where we find no evidence for it. 

A further eleven galaxies show DIP values greater than 0.5 at the 1 to 2 a level — i.e., to have 
potential deviations from unimodality. These are NGC 584, NGC 596, NGC 821, NGC 1023, 
NGC 2434, NGC 2778, NGC 3115, NGC 3610, NGC 4365, NGC 4458, and NGC 4478. The studies 
of Neilsen k Tsvetanov (1999) and Kissler-Patig et al. (1999) show that NGC 4478 does not seem 
to have two clearly distinct populations. In contrast, the recent study of Kundu & Whitmore 
(1998) reveals that NGC 3115 hosts two distinct populations. We speculate that at least some of 
these systems might have several distinct populations. Further detailed analysis are required for 
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these candidates. 

The remaining 23 galaxies do not show deviations from unimodality. We note that the 
detection of a dip does not correlate with the number of globular clusters detected; the sample 
size is taken into account when computing the significance of the statistic. 

In Figure 2, the bimodality plots of interest are those which show DIP significance versus 
both the absolute magnitude and the environment density of the host galaxy. The detection 
of bimodality does not seem to correlate with the environment. However, some trend can be 
seen with absolute magnitude: low-luminosity systems (Mv > —20) appear unimodal. This 
trend has already been noted by Kissler-Patig (1997) in relation to other differences between 
low- and high-luminosity systems. It is not clear yet how much of these differences can be 
associated with fundamental differences between the two types of galaxies, and how much is 
just a bias in the observations. Indeed, if the blue peak is constant in all galaxies, while the 
red one scales with galaxy luminosity, the two peaks will be harder to disentangle in fainter 
galaxies (cf. Kissler-Patig et al. 1998a). On the other hand, high-luminosity systems do not 
appear systematically to show bimodality. We note the interesting cases of NGC 4406 (M86) and 
NGC 4374, two giant ellipticals in Virgo. The first galaxy, a central giant elliptical dominating 
an infalling group (Schindler et al. 1999), is clearly unimodal. NGC 4374 also has a unimodal 
distribution; however, it is heavily skewed to the red, so that a large number of red (presumably 
metal-rich) globular clusters exist but seem to form a continuous distribution rather than a 
distinct population. 

In summary, it seems harder to rule out the presence of several populations than to detect 
them. In our sample, roughly half of the systems show multiple populations of globular clusters, 
independently of the environment but preferentially in the brightest galaxies. Our future challenge 
will be to confirm or rule out single populations in low-luminosity systems. 



7. Combining Galaxy Samples 

The moments of the V — I distributions are not well estimated in many of our galaxies 
based on small numbers of globular clusters. To compare correlations more robustly, we combine 
the different globular cluster samples to give better estimates of the moments. Because we 
do not have to exclude galaxies based on there being a small number of measured globular 
clusters, we now have 43 galaxies in the combined sample. We use the four main parent galaxy 
properties — magnitude, a, Mg2 index, and density environment — for this analysis, breaking the 
sample into two bins. Figure 5 plots the results for the combined V — I distributions. Table 4 
presents the various cuts used for the combined samples and the moments and DIP probability 
for each sub-sample. Table 4 also presents the number of galaxies and the number of globular 
clusters in each sub-sample. Each sub-sample has greater than 800 clusters where the noise 
from the moment estimation is minimal. The choice for the cuts in the sub-samples results from 
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our inspection of the particular distribution; we attempted to find a clear division between the 
samples, and, when none was apparent, we used cuts which had roughly equal range among the 
particular variable. The cuts are partly subjective, but the results do not depend strongly on 
them. 

The central location of the V — I peak is statistically different only in the division between 
faint and bright galaxies. The faint galaxies have a peak around 0.94, and the bright galaxies have 
a peak around 1.13. For the other divisions, the V — I peaks are similar, but the tail indexes differ 
significantly. Either faint galaxies, low a, low Mg2, or low density cause the V — I distribution to 
have significantly more tail weight than compared to their counterpart sample. 

The most extreme differences between the combined samples results from the DIP probability. 
Galaxies that are faint, have low a or low Mg2 are consistent with a unimodal distribution, 
whereas their counterpart samples are highly non-unimodal. This effect is also seen when looking 
at individual galaxies (see Section 6 above), and it appears to be a fundamental property of 
galaxies. 



8. Discussion 

8.1. The Color Distributions 

An important discovery in the study of globular clusters systems in the last decade is the 
fact that the globular clusters' color distribution of many ellipticals show several peaks (Zepf & 
Ashman 1993). This property has been interpreted as the presence of several distinct populations 
of globular clusters around the host galaxy. The origin of these different populations is still under 
debate and might well have several causes. The original explanation argues that spiral-spiral 
mergers provide old globular clusters from the progenitors as well as produce a new population 
during the merger (Ashman & Zepf 1992, see also Schweizer 1987). However, the accretion of 
enough dwarf galaxies with metal poor clusters will also obviously build up a distinct metal-poor 
population. Cote et al. (1998), and Hilker (1999) investigated this idea. Yet another process 
involves the pre-galactic formation of metal-poor clusters as the building blocks of galaxies, 
followed by the formation of more metal-rich globular clusters during the collapse of the galaxy 
(see Kissler-Patig 1997, and Kissler-Patig et al. 1998a). Only slightly different is the theory that 
galaxies formed all their globular clusters in situ and had two main phases of star formation 
(e.g. Forbes et al. 1997, Harris et al. 1998, Harris et al. 1999). Today's challenge is to discriminate 
between the various scenarios with detailed comparison to the present data. 

The diversity of the color distributions present in our sample indicates that indeed the 
scenario is complex. We discuss in turn some of the groups: 

• Galaxies that host "only red clusters" — i.e., those having a red peak and a small scale value. 
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The existence of a blue population in these galaxies is not excluded but is unlikely to be significant. 
The favored scenario for these galaxies is that a single dominant collapse created a bulge and a 
red population of clusters. We cannot rule the accretion scenario out if these galaxies were formed 
in an environment poor in material suitable for accretion (e.g. dwarf galaxies, cf. Chapelon 1999). 
These galaxies would then be in low density environments or in the outskirts of clusters; neither is 
confirmed nor ruled out by our present statistics. Finally, the merger scenario could also work for 
these galaxies if the progenitors were poor in blue clusters. If the merger happened at early times, 
the progenitors could be still in a gaseous form, or, alternatively, could be low-surface-brightness 
galaxies for late mergers in the field that host a large amount of gas but apparently few clusters. 
However, a scenario with gaseous mergers at early times is hard to distinguish from the two-stage 
in situ scenario. 



• Galaxies that host "only blue clusters" — i.e., that have a blue peak and a small scale 
value. These galaxies do not seem to host a significant population of red (presumably metal-rich) 
clusters. In this case, none of the above scenarios would easily be able to explain such systems; one 
would have to advocate the formation of bulge stars without the formation of clusters. More likely 
is a "blue" red peak, either because the red clusters formed with a low metallicity or because they 
have a significant younger age (about half as old as the true blue clusters, cf. Kissler-Patig et al. 
1998a) that make them appear blue. In these last two cases, all scenarios work, with some 
preference to the merger scenario in the frame of hierarchical clustering models that predict late 
mergers in the field. 

• Galaxies with broad or bimodal color distributions — i.e., that which presumably host a 
comparable amount (say within a factor of about 3) of blue and red clusters. These are the cases 
for which the scenarios described above were designed, and we refer to the respective papers for 
the pros and cons. 

Finally, we note that as long as some information is missing (merger epoch, globular clusters 
formation mechanisms, faint end of the luminosity function in different environments, etc.), and 
given the diversity of scenarios, it seems that not always the same mechanism is the dominant one. 



8.2. Field and Cluster Systems 

Although based on small number statistics, field and cluster early-type galaxies seem to differ 
in the properties of their globular cluster systems. Whereas cluster galaxies show a correlation 
between peak value and the galaxy size (as measured by the velocity dispersion), such a relation 
must be ruled out for field galaxies. 

For cluster ellipticals, the recent spectroscopic studies (Cohen et al. 1998, Kissler-Patig et al. 
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1998b) show that all globular clusters are old (around 10 Gyr or older). At these ages, the clusters 
can be considered coeval with respect to the influence of their ages on color. Color differences can 
therefore be interpreted as almost pure metallicity differences, and the relation between peak color 
and galaxy luminosity might reflect the equivalent of the Mg2-a relation for galaxies. Forbes et al. 
(1997) suggest that in systems with two distinct globular cluster populations, the mean metallicity 
of the red clusters correlates with the luminosity of the host galaxy; the blue peak does not. Their 
conclusion implies that the correlation of the total sample's mean color with galaxy size is driven 
by the metal-rich globular clusters only. Consequently, this would support scenarios in which the 
metal-poor globular clusters formed in events uncorrelated with their future host galaxy — e.g., in 
systems that are separate entities from their future hosts. On the other hand, the same research 
suggests that the metal-rich globular clusters formed in events correlated with the future host 
galaxy — e.g., the formation of a bulge in a collapse or merger event. In addition, this correlation 
suggests that the blue sub-populations have roughly a constant color, while for the red population 
both the peak shifts and the scale increase as a function of color. 

Perhaps surprisingly, the peak color of globular clusters in field galaxies does not correlate 
with the velocity dispersion (~ size) of their host galaxy. This result appears to be in contrast with 
the observation that field and cluster galaxies do exhibit the same Mg2-o" relations (Bernardi et al. 
1998). We see several possible explanations. 

A first alternative might be that the correlation for cluster galaxies is artificial. We verified 
that it is not caused by sampling different regions in the different galaxies. Indeed, the ratio of blue 
to red galaxies does vary with radius (Geisler et al. 1997, Kissler-Patig et al. 1997b, Puzia et al. 
1999). Furthermore, the WFPC2 field samples the outer parts less in larger galaxies than in small 
ones, so that the relative ratio of blue to red galaxies is systematically smaller in large galaxies at 
fixed field of view. Puzia et al. (1999) tested the importance of this effect on NGC 4472, which 
has exposures centered on the galaxy as well as several arcminute offsets north and south. The 
respective peaks of the blue and red populations do not change with radius within the errors; 
however, the varying ratio of blue to red causes the overall peak to shift by 0.02 mag from the 
inner to the outer field. The effect is too small, though, to explain any systematic correlation. 
Moreover, in our sample, field and cluster galaxies span different values of galaxy a. Where cluster 
galaxies are of intermediate to giant size, the field galaxies span the low to intermediate size range. 
If indeed intermediate and giant early-type galaxies differ (see above, cf. Kissler-Patig 1997), the 
presence/absence of a correlation would be due not to environment, but rather to a function of the 
galaxy size. For example, one could speculate that small galaxies are not surrounded by a "halo" 
of blue clusters, thus showing a small trend of system peak with size. Intermediate size or larger 
galaxies might or might not have such a halo of blue clusters, complicating the interpretation 
of the peak color and introducing a scatter due to the possible existence of blue clusters. Giant 
ellipticals would always have blue clusters, and a significant redder population that correlates 
with galaxy size, leading to an overall trend of peak color with galaxy size again. In other words, 
the correlation would be an artifact of the mix of populations which show a correlation with 
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environment, rather than a physical difference in the enrichment histories. 

A second alternative would be a physical interpretation of the result. Thus, the presence of 
a Mg2-cj relations in field and cluster galaxies, but absence/presence respectively of a peak-o" 
relation, would indicate that stars and globular clusters are coupled in cluster galaxies but not in 
field galaxies. We regard this alternative as unlikely, however, since cluster formation and star 
formation appear to go together. 

Finally, Mg2 and V — I color could vary differently with age and metallicity. The variations 
of Mg2 with age and metallicity are still uncertain and if, for example, the metal-rich population 
in field ellipticals show a large spread in age from galaxy to galaxy, a trend could mimic an 
Mg2-(j relation in the field but not a peak-a relation. In contrast, since all globular clusters in 
cluster galaxies are old, such a trend would not be present, and both Mg2 and peak would trace 
metallicity only. Clearly such speculations require a much better understanding of the influence 
of age and metallicity on the properties of composite stellar populations. Detailed studies of the 
ages and metallicities of the systems are needed. 



9. Summary and Outlook 

The results presented in this paper suggest that we can use globular cluster systems to study 
the merger/formation history of both itself and its host galaxy. The differences we see between 
the field and cluster systems suggest that the cluster galaxies are more influenced by accretion, 
whereas the field systems may be the result of "in situ" globular cluster formation. Combining 
this result with the strong signature of multi-modality seen in many of the systems provides an 
observational comparison to the various formation scenarios. In a companion paper (Cote et al. 
2000), we will compare these results to the theoretical predictions of Cote et al. (1998). More 
detailed studies of the systems that are candidates for bimodal distributions are underway, 
including follow-up near-infrared photometry and multi-object spectroscopy. 

All photometry results are available from the authors. 

We are indebted to Pat Cote for his help and advises in some of the early-stages of this 
study. Thanks to Thomas Puzia for his help with the NGC 4472 data and extra-tests. We thank 
Alvio Renzini for enlightening discussions, as well as the ENEAR team for providing galaxy data 
prior to publication, especially Mariangela Bernardi for her help in digging these latter out. KG 
is supported by NASA through Hubble Fellowship grant HF-01090.01-97A awarded by the Space 
Telescope Science Institute, which is operated by the Association of the Universities for Research 
in Astronomy, Inc., for NASA under contract NAS 5-26555. MKP was partly supported by a 
Feodor Lynen Fellowship of the Alexander von Humboldt Foundation. 
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Table 1. Observational Parameters 



Galaxy 


Filter 


T 

J- cxp 


N 


PI (number) 


N0584 


F555W 


1520 


5 


Faber (6099) 




F814W 


1450 


7 


Faber (6099) 


N0596 


F555W 


1600 


4 


Richstone (6587) 




F814W 


1800 


3 


Richstone (6587) 


N0821 


F555W 


1520 


5 


Faber (6099) 




F814W 


1450 


7 


Faber (6099) 


N1023^ 


F555W 


2920 


11 


Faber.Richstonc (6099,6587) 




F814W 


3780 


10 


Faber,Riclistone (6099,6587) 


N1316 


F450W 


5000 


5 


GriUmair (5990) 




F814W 


1860 


3 


GriUmair (5990) 


N1399 


F450W 


5200 


4 


GriUmair (5990) 




F814W 


1800 


3 


GriUmair (5990) 


N1404 


F450W 


5000 


5 


GriUmair (5990) 




F814W 


1860 


3 


GriUmair (5990) 


N1426 


F555W 


1600 


4 


Richstone (6587) 




F814W 


1800 


3 


Richstone (6587) 


N1700 


F555W 


1600 


4 


Richstone (6587) 




F814W 


1800 


3 


Richstone (6587) 


N2300 


F555W 


1520 


5 


Faber (6099) 




F814W 


1450 


7 


Faber (6099) 


N2434 


F555W 


1300 


2 


CaroUo (5943) 




F814W 


1000 


2 


CaroUo (5943) 


N2778 


F555W 


1520 


5 


Faber (6099) 




F814W 


1450 


7 


Faber (6099) 


N3115 


F555W 


1015 


6 


Faber (5512) 




F814W 


1236 


6 


Faber (5512) 


N3377 


F555W 


1330 


6 


Faber (5512) 




F814W 


1236 


6 


Faber (5512) 


N3379 


F555W 


1660 


4 


Faber (5512) 




F814W 


1340 


4 


Faber (5512) 


N3384 


F555W 


1430 


6 


Faber (5512) 




F814W 


1290 


6 


Faber (5512) 


N3585 


F555W 


1200 


6 


Richstone (6587) 




F814W 


1900 


3 


Richstone (6587) 


N3610 


F555W 


1400 


6 


Richstone (6587) 




F814W 


2100 


3 


Richstone (6587) 


N3640 


F555W 


1200 


6 


Richstone (6587) 




F814W 


2100 


3 


Richstone (6587) 


N4125 


F555W 


1400 


6 


Faber (6587) 




F814W 


2100 


3 


Faber (6587) 


N4192 


F555W 


660 


3 


Rubin (5375) 




F814W 


660 


3 


Rubin (5375) 


N4291 


F555W 


1520 


5 


Faber (6099) 




F814W 


1450 


7 


Faber (6099) 


N4343 


F555W 


660 


3 


Rubin (5375) 
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Table 1 — Continued 



Galaxy 


Filter 


J- cxp 


Noxp 


PI (number) 




F814W 


660 


3 


Rubin (5375) 


N4365 


F555W 


2200 


2 


Brodie (5920) 




F814W 


2300 


2 


Brodie (5920) 


N4374 


F547M 


1200 


2 


Bower (6094) 




F814W 


520 


2 


Bower (6094) 


N4406 


F555W 


1500 


3 


Faber (5512) 




F814W 


1500 


3 


Faber (5512) 


N4450 


F555W 


520 


3 


Rubin (5375) 




F814W 


520 


3 


Rubin (5375) 


N4458 


F555W 


1340 


4 


Faber (5512) 




F814W 


1120 


5 


Faber (5512) 


N4472'' 


F555W 


4400 


4 


Brodie (5920) 




F814W 


4600 


4 


Brodie (5920) 


N4473 


F555W 


1800 


3 


Faber (6099) 




F814W 


2000 


4 


Faber (6099) 


N4478 


F555W 


1600 


4 


Richstone (6587) 




F814W 


1800 


3 


Richstone (6587) 


N4486 


F555W 


2810 


6 


Macchetto (5477) 




F814W 


2430 


5 


Macchetto (5477) 


N4486B 


F555W 


1800 


3 


Faber (6099) 




F814W 


2000 


4 


Faber (6099) 


N4526 


F555W 


520 


3 


Rubin (5375) 




F814W 


520 


3 


Rubin (5375) 


N4536 


F555W 


660 


3 


Rubin (5375) 




F814W 


660 


3 


Rubin (5375) 


N4550 


F555W 


1200 


3 


Rubin (5375) 




F814W 


1200 


3 


Rubin (5375) 


N4552 


F555W 


2400 


4 


Faber (6099) 




F814W 


1500 


3 


Faber (6099) 


N4569 


F555W 


520 


3 


Rubin (5375) 




F814W 


520 


3 


Rubin (5375) 


N4594 


F547M 


1340 


4 


Faber (5512) 




F814W 


1470 


6 


Faber (5512) 


N4621 


F555W 


1380 


6 


Faber (5512) 




F814W 


1280 


6 


Faber (5512) 


N4649 


F555W 


2100 


2 


Westphal (6286) 




F814W 


2500 


2 


Westphal (6286) 


N4660 


F555W 


1000 


5 


Faber (5512) 




F814W 


850 


5 


Faber (5512) 


N4881 


F555W 


7200 


8 


Westphal (5233) 




F814W 


7200 


8 


Westphal (5233) 


N5018 


F555W 


1200 


6 


Richstone (6587) 




F814W 


1900 


3 


Richstone (6587) 


N5061 


F555W 


1200 


6 


Richstone (6587) 




F814W 


1900 


3 


Richstone (6587) 
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Table 1 — Continued 



Galaxy 


Filter 


J- cxp 


Nexp 


PI (number) 


N5845 


F555W 


2140 


5 


Fabcr (6099) 




F814W 


1120 


5 


Fabor (6099) 


N5846'' 


F555W 


6600 


6 


Brodie (5920) 




F814W 


6900 


6 


Brodie (5920) 


N7192 


F555W 


1300 


2 


CaroUo (5943) 




F814W 


1000 


2 


CaroUo (5943) 


N7457 


F555W 


1380 


6 


Faber (5512) 




F814W 


1280 


6 


Faber (5512) 


IC4889 


F555W 


600 


1 


CaroUo (5943) 




F814W 


600 


1 


CaroUo (5943) 



''N1023 and N4472 have two different pointings; 
N5846 has three. 
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Table 2. Galaxy Properties 



Galaxy 


Type 


Mv 


log(o-) 


Mg2 


log(p) 


a4 


log(Re)(pc) 


N584 


E4 


-21.16 


2.34 


0.29 


0.42 


1.50 


3.51 


N596 


Epec 


-20.93 


2.22 


0.25 


0.40 


1.30 




N821 


E6 


-20.84 


2.32 


0.32 


0.08 


2.50 


3.45 


N1023 


SBO 


-20.61 


2.34 


0.34 


0.57 




3.26 


N1399 


Elpec 


-21.97 


2.52 


0.37 


1.59 


0.00 




N1404 


El 


-21.56 


2.38 


0.34 


1.59 


0.50 




N1426 


E4 


-20.42 


2.18 


0.28 


0.66 




3.44 


N1700 


E4 


-21.90 


2.36 


0.28 




0.90 




N2300 


SAO 


-20.76 


2.42 


0.33 


0.14 


-0.60 


3.58 


N2434 


EO 


-20.60 


2.31 


0.27 


0.19 




3.49 


N2778 


E 


-19.43 


2.26 


0.34 


0.28 






N3115 


SO 


-20.82 


2.45 


0.34 


0.08 


8.00 


3.17 


N3377 


E5 


-19.96 


2.18 


0.29 


0.49 


1.20 


3.21 


N3379 


El 


-20.83 


2.35 


0.34 


0.52 


0.10 


3.23 


N3384 


SBO 


-19.89 


2.21 


0.29 


0.54 






N3585 


E7 


-21.75 


2.34 


0.32 


0.12 


4.50 


3.50 


N3610 


E5 


-20.92 


2.21 


0.27 


0.30 


2.50 




N3640 


E3 


-21.75 


2.24 


0.27 


0.18 


-0.20 




N4125 


E6pec 


-22.09 


2.25 


0.32 


0.34 


1.00 


3.43 


N4291 


E 


-20.69 


2.43 


0.32 


0.36 


-0.40 


3.64 


N4365 


E3 


-22.12 


2.42 


0.35 


2.93 


-1.10 


3.79 


N4374 


El 


-22.21 


2.44 


0.32 


3.99 


-0.40 




N4406 


SO 


-22.41 


2.40 


0.33 


1.41 


-0.70 




N4458 


EO 


-19.06 


2.02 


0.25 


3.21 






N4472 


E2 


-22.59 


2.48 


0.34 


3.31 


-0.30 


3.89 


N4473 


E5 


-20.80 


2.26 


0.33 


2.17 


0.90 


3.60 


N4478 


E2 


-19.94 


2.13 


0.27 


3.92 


-0.80 




N4486 


EO 


-22.48 


2.56 


0.30 


4.17 


0.00 


3.89 


N4486b 


cEO 


-17.75 


2.30 


0.30 


4.17 




3.51 


N4526 


SABO 


-20.89 


2.47 


0.30 


2.45 




3.57 


N4550 


SBO 


-18.81 


2.01 


0.17 


2.97 


-0.70 




N4552 


E 


-21.26 


2.42 


0.35 


2.97 




3.35 


N4621 


E5 


-21.63 


2.40 


0.35 


2.60 


1.50 


3.54 


N4649 


E2 


-22.31 


2.56 


0.37 


3.49 


-0.50 


3.74 


N4660 


E 


-19.26 


2.29 


0.33 


3.37 


2.70 


3.56 


N4881 


E 


-21.60 


2.33 


0.29 








N5018 


E3 


-22.31 


2.32 


0.21 


0.29 


2.00 




N5061 


EO 


-21.62 


2.30 


0.27 


0.31 


-1.80 


3.90 


N5845 


E 


-19.63 


2.42 


0.32 


0.84 


0.80 




N5846 


EO 


-22.23 


2.40 


0.34 


0.84 


0.00 


3.37 


N7192 


SAO 


-21.49 


2.27 


0.26 


0.28 




3.71 


N7457 


SAO 


-18.25 


1.76 


0.19 


0.13 






IC4889 


E5 


-21.26 


2.22 


0.26 


0.13 


1.00 
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Table 3. Moments oi V — I Distributions 



Galaxy Ngc V — I peak 



Scale 



Skewness 



Tail 



DIP 



N584 


67 


1.092 


N596 


31 


1.049 


N821 


55 


1.008 


N1023 


149 


1.037 


N1316 


186 


0.966 


N1399 


410 


1.051 


N1404 


138 


1.017 


N1426 


74 


1.084 


N1700 


27 


1.122 


N2300 


82 


1.071 


N2434 


143 


0.945 


N2778 


18 


1.072 


N3115 


94 


1.030 


N3377 


86 


1.072 


N3379 


52 


1.095 


N3384 


39 


1.023 


N3585 


70 


1.110 


N3610 


29 


1.105 


N3640 


48 


1.080 


N4125 


68 


1.001 


Nil 92 


49 


0.884 


N4291 


71 


0.983 


N4343 


25 


1.026 


N4365 


409 


1.094 


N4374 


123 


0.999 


N4406 


131 


1.041 


N445O 


29 


0.941 


N4458 


33 


1.071 


N4472 


486 


1.095 


N4473 


112 


1.026 


N4478 


50 


1.060 


N4486 


655 


1.108 


N4486b 


87 


0.978 


N4526 


75 


1.030 


N4536 


87 


0.766 


N4550 


37 


1.020 


N4552 


171 


1.062 


N4569 


67 


0.816 


N4594 


122 


1.065 


N4621 


126 


1.079 


N4649 


378 


1.096 


N4660 


74 


0.962 


N4881 


84 


0.915 


N5018 


10 


0.931 


N5061 


84 


1.099 



1.123 
1.062 
1.096 
1.004 
1.041 
0.976 
1.063 
1.010 
0.981 
0.951 
1.063 
1.040 
1.037 
0.999 
1.118 
1.050 
1.187 
1.059 
1.110 
1.031 
0.968 
0.922 
1.143 
1.007 
1.058 
1.002 
1.095 
1.047 
1.129 
1.054 
1.084 
0.969 
1.149 
1.071 
1.182 
1.044 
1.116 
1.045 
1.042 
0.954 
0.970 
0.804 
1.016 
0.952 
1.096 
0.962 
1.111 
1.079 
1.021 
0.977 
1.059 
1.021 
1.017 
0.850 
1.139 
0.982 
1.113 
1.078 
1.050 
1.003 
1.100 
1.022 
1.118 
1.097 
1.001 
0.953 
1.074 
0.988 
0.847 
0.691 
1.098 
0.944 
1.083 
1.043 
0.885 
0.748 
1.089 
1.042 
1.099 
1.059 
1.111 
1.081 
0.992 
0.935 
0.959 
0.868 
1.013 
0.848 
1.127 
1.068 



0.159 
0.158 
0.142 
0.198 
0.130 
0.124 
0.137 
0.176 
0.213 
0.193 
0.189 
0.213 
0.167 
0.141 
0.172 
0.204 
0.202 
0.191 
0.158 
0.218 
0.366 
0.161 
0.246 
0.188 
0.151 
0.132 
0.284 
0.221 
0.250 
0.149 
0.166 
0.158 
0.142 
0.214 
0.388 
0.305 
0.156 
0.347 
0.152 
0.128 
0.185 
0.148 
0.276 
0.187 
0.190 



0.192 
0.130 
0.201 
0.115 
0.158 
0.126 
0.229 
0.172 
0.149 
0.114 
0.130 
0.118 
0.152 
0.125 
0.198 
0.153 
0.292 
0.163 
0.222 
0.168 
0.224 
0.156 
0.306 
0.138 
0.190 
0.146 
0.169 
0.118 
0.208 
0.142 
0.261 
0.157 
0.239 
0.173 
0.245 
0.149 
0.186 
0.132 
0.251 
0.185 
0.425 
0.306 
0.187 
0.138 
0.325 
0.137 
0.203 
0.175 
0.171 
0.133 
0.146 
0.119 
0.373 
0.189 
0.270 
0.127 
0.269 
0.232 
0.168 
0.134 
0.206 
0.136 
0.164 
0.151 
0.175 
0.117 
0.248 
0.186 
0.439 
0.332 
0.363 
0.254 
0.171 
0.142 
0.390 
0.309 
0.176 
0.130 
0.141 
0.116 
0.200 
0.171 
0.183 
0.117 
0.337 
0.225 
0.309 
0.111 
0.242 
0.143 



-0.013 
+0.306 
+0.257 
+1.660 
-1.001 
-4.387 
+0.125 
+0.486 
+1.328 
+0.881 
+0.465 
-0.827 
+1.363 
-0.680 
+0.896 
+1.248 
+0.331 
+0.154 
+0.982 
+0.241 
+0.179 
+0.637 
+0.586 
+0.253 
+2.025 
+0.283 
-0.278 
+0.887 
-0.236 
-0.618 
+0.267 
+0.290 
+0.501 
-0.243 
+2.334 
+0.022 
+ 1.121 
+0.492 
-1.988 
-1.611 
-1.401 
+0.899 
+0.505 
+0.193 
-0.473 



+ 1.128 
-1.159 
+ 1.541 
-0.981 
+ 1.518 
-1.006 
+2.953 
+0.271 
+0.499 
-2.382 
-2.187 
-6.649 
+ 1.741 
-1.470 
+ 1.660 
-0.673 
+2.416 
+0.034 
+2.135 
-0.377 
+ 1.706 
-0.723 
+0.597 
-2.033 
+2.652 
+0.119 
+0.523 
-1.905 
+2.078 
-0.307 
+2.285 
+0.058 
+ 1.564 
-0.971 
+ 1.467 
-1.102 
+2.125 
-0.250 
+ 1.560 
-1.019 
+ 1.268 
-0.910 
+ 1.770 
-0.561 
+ 1.667 
-0.622 
+ 1.851 
-1.335 
+3.285 
+0.560 
+ 1.730 
-1.347 
+0.889 
-1.476 
+ 1.970 
-0.380 
+ 1.308 
-1.865 
+0.983 
-1.992 
+ 1.492 
-1.054 
+2.615 
-2.404 
+ 1.678 
-0.746 
+0.993 
-1.480 
+3.478 
+0.951 
+ 1.218 
-1.115 
+2.603 
-0.354 
+ 1.667 
-0.798 
-0.430 
-3.310 
-0.130 
-3.073 
+0.251 
-2.909 
+2.148 
-0.330 
+ 1.683 
-0.497 
+ 1.340 
-1.145 
+0.723 
-1.602 



1.080 
1.176 
0.805 
1.007 
1.196 
0.886 
0.770 
0.965 
1.093 
0.911 
1.100 
1.291 
0.933 
1.017 
0.950 
1.067 
0.993 
0.956 
0.852 
1.007 
0.888 
1.072 
1.442 
0.941 
0.909 
0.777 
1.360 
0.928 
0.985 
0.896 
1.017 
0.821 
1.038 
0.815 
1.042 
0.909 
0.959 
0.988 
0.848 
0.916 
0.867 
1.091 
1.321 
1.714 
1.231 



1.484 
0.864 
2.049 
0.833 
0.980 
0.690 
1.401 
0.878 
1.479 
0.971 
0.998 
0.784 
0.878 
0.693 
1.205 
0.817 
1.937 
0.778 
1.110 
0.772 
1.473 
0.875 
2.567 
0.772 
1.119 
0.773 
1.250 
0.826 
1.317 
0.728 
1.640 
0.746 
1.276 
0.805 
1.513 
0.747 
1.283 
0.697 
1.378 
0.784 
1.254 
0.680 
1.348 
0.865 
2.543 
0.863 
1.047 
0.859 
1.111 
0.773 
0.923 
0.690 
2.018 
0.867 
1.937 
0.556 
1.087 
0.905 
1.052 
0.765 
1.402 
0.797 
0.881 
0.764 
1.491 
0.844 
1.037 
0.644 
1.329 
0.836 
1.238 
0.700 
1.094 
0.851 
1.181 
0.807 
1.121 
0.687 
1.119 
0.773 
0.946 
0.794 
1.445 
0.832 
1.770 
1.026 
3.228 
0.965 
1.812 
0.907 



0.649 
0.661 
0.785 
0.635 
0.500 
0.968 
0.995 
0.995 
0.500 
0.500 
0.598 
0.595 
0.709 
0.500 
0.500 
0.500 
0.500 
0.996 
0.971 
0.500 
0.500 
0.500 
0.500 
0.757 
0.520 
0.500 
0.500 
0.698 
0.997 
0.500 
0.614 
0.978 
0.500 
0.999 
0.500 
0.500 
0.500 
0.973 
0.995 
0.539 
0.961 
0.500 
0.500 
0.500 
0.500 



0.852 
0.500 
0.778 
0.500 
0.840 
0.500 
0.696 
0.500 
0.590 
0.500 
0.969 
0.908 
0.995 
0.886 
1.000 
0.705 
0.743 
0.500 
0.675 
0.500 
0.798 
0.500 
0.800 
0.500 
0.754 
0.500 
0.779 
0.500 
0.770 
0.500 
0.737 
0.500 
0.768 
0.500 
1.000 
0.557 
0.994 
0.593 
0.814 
0.500 
0.668 
0.500 
0.941 
0.500 
0.827 
0.500 
0.795 
0.514 
0.690 
0.500 
0.571 
0.500 
0.768 
0.500 
0.764 
0.500 
0.997 
0.945 
0.667 
0.500 
0.744 
0.500 
0.978 
0.951 
0.767 
0.500 
0.999 
0.500 
0.840 
0.500 
0.688 
0.500 
0.707 
0.500 
0.988 
0.554 
1.000 
0.627 
0.654 
0.500 
0.963 
0.810 
0.767 
0.500 
0.999 
0.500 
0.773 
0.500 
0.800 
0.500 
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Table 3 — Continued 



Galaxy Ngc V — I peak 



Scale 



Skewness 



Tail 



DIP 



N5845 
N5846 
N7192 
N7457 
IC4889 



28 
781 
95 
48 
48 



1.127 lill 
1.080 ITe 
1.068 nil 
1.024 lilt 
1.024 i-SIJi 



0.364 
0.211 
0.266 
0.232 
0.313 
0.215 

0.147 g:??i 



0.285 
0.249 
0.262 



0.147 



0.203 
0.113 



+0.543 
-2.895 
-0.196 
+0.673 
+0.673 



+ 1.569 
-0.617 
-1.318 
-4.551 
+0.990 
-1.338 
+ 1.802 
-0.461 
+ 1.800 
-0.530 



1 070 1.851 

l-^'^ 0.878 

1-J^O/ 1070 

1 1 1 p- 1.446 

J^-J^J^O 0.882 

1 oae 2.449 

i.28» 0.857 

1 OQQ 2.521 

-l-^oo 0.875 



0.500 gilgg 
0.895 °;ii 
0.500 giUg 
0.500 gii^g 

0.500 g:igg 



Table 4. Moments of Combined Samples 



Sub-Sample 



Range 



Nr 



Ngc V — I peak 



Scale 



Skewness 



Tail 



DIP 



Faint Mv > -20.70 15 803 

Bright M^. < -21.75 13 2737 



+ 1 663 +^°°^ 

-1-1. UOO ^0.302 



-1.397 



+0.518 
3.383 



0.941 

0.884 Ulr 



0.954 "oTai 



0.500 
0.963 



Low a 
High a 



log((7) < 2.3 
log(o-) > 2.4 



16 
12 



872 
2992 



1.041 nil 

1.075 ;:g?g 



0.188 r,% 

r, 170 0.177 
U.4' J 0.169 



+ 1.243 ±^:?« 
-0.455 li-fj 



1 r|99 1.080 
i.UZZ 0.964 

0.879 g;ig^ 



0.759 °o:fot 
0.996 g:^^ 



Low Mg2 
High Mg2 



Mg2 < 0.30 
Mg2 > 0.34 



18 
15 



1496 
2597 



1.069 
1.070 



0.205 o!i97 
0.170 nil 



-0.104 
-0.292 



+ 1.786 
-1.894 
+ 1.618 
-2.258 



1.039 S:^^^ 



0.791 



0.888 gJs^ 0.985 



Low p 
High p 



log(p) < 1.0 
log(p) > 1.0 



24 
17 



2202 
3441 



1.057 ntt 
1.069 lilt 



0.200 "oigl 
0.170 °oil6 



+1.306 
-0.262 



0.967 : 



0.887 !]:i|31 0.999 
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Fig. 1. — Histogram of the globular cluster system's V — I colors for all of the galaxies. The solid 
line represents a density estimate and the dotted lines are the 90% confidence bands. The total 
number of globular clusters appears under the galaxy name. 
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Fig. 1. — continued. 
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Fig. 1. — continued. 
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Fig. 1. — continued. 
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Fig. 2. — Plots of various globular cluster systems properties and parent galaxy properties. The 
triangles represent field galaxies and the circles represent cluster galaxies. The names along the 
diagonal are the variable which is plotted along that particular row and column. Each plot is 
repeated twice here as the plots in the lower-left triangle and the upper-right triangle represent the 
same points with the axes inverted. 
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Fig. 3. — Scale and log(cr) versus the V — I color of the globular cluster systems. The colors of the 
points represent the same as in the previous plot. 
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Fig. 4. — The distribution of the DIP probabiHty with it's 
galaxies. The DIP probabihty is bounded by 0.5 and 1.0. 
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Fig. 5. — V — I distributions and their 90% confidence bands for the globular clusters combined 
across galaxies; we use two bins for each of the named galaxy properties. Each curve represents an 
adaptive kernel estimate of the density with unity area. 



